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SUMMARY 

The sky subtraction performances of multi-fiber spectrographs are discussed, 
analysing in detail the case of the OPTOPUS system at the 3.6 meter ESO 
telescope at La Silla. A standard technique, based on flat-fields obtained with 
a uniformly illuminated screen on the dome, provides poor results. A new method 
has been developed, using the [01] emission line at 5577 A as a calibrator of the 
fiber transmittance, taking into account the diffuse light and the influence of each 
fiber on the adjacent ones, and correcting for the effects of the image distortions 
on the sky sampling. In this way the accuracy of the sky subtraction improves 
from 2 - 8% to 1.3 - 1.6%. 

1. INTRODUCTION 

Spectrographs fed by optical fibers allow to observe many (from tens to 
hundreds) objects simultaneously. The higher the spatial density of the targets, 
the larger the advantage over "classic" single-slit spectroscopy, thus a considerable 
effort has been spent in extending fiber spectroscopy to fainter and fainter objects. 
To this end and in view of the application of this technique to telescopes of 
the eight-ten meter class, the limitations imposed by the accuracy in the sky 
subtraction have become a key issue. 

Various papers (Elston and Barden 1989, Wyse and Gilmore 1992, Mignoli 
and Cuby 1994) have discussed this point with the following results: 

1) An accuracy of the sky subtraction in the range 1 — 2% is considered good 
in practice, while, on theoretical grounds, a limiting accuracy of a percent or 
better is expected to be achievable. 

2) Such a limitation is due to: focal ratio degradation of the fibers, internal 
scattered light, sky variations across the field (caused by blue galaxies, airglow, 
etc.), telecentricity t effects (Wynne 1993), contamination from adjacent 
fibers (cross-talk), poor determination of the fiber transmittance and of the 
wavelength calibration. 

t Telecentricity: effect due to the misalignment between the axis of the optical 
fibers and the light beam axis coming from the telescope output pupil. 
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Section 2 describes a "standard" approach to the sky subtraction and how to 
measure the accuracy of the method. In section 3 the importance of the various 
sources of noise is analyzed. Section 4 outhnes the basics of a new sky subtraction 
method and presents the comparison with the results of section 3. Finally, section 
5 deals with further problems regarding the extraction of the spectra. 

2. ESTABLISHING THE ACCURACY OF THE SKY SUBTRAC- 
TION 

2.1) the instrument 

In the following we will analyze the sky subtraction performances of a 
multifiber spectrograph in the case of the OPTOPUS system operating at the 
3.6 meter telescope of ESQ La Silla, described by Avila and D'Odorico (1993) and 
in the ESQ Users Manual (Schwarz and Melnick, 1993). The instrument consists 
of 50 fibers, manually plugged in a previously drilled metal template and feeding a 
B. & C. grating spectrograph. Each fiber has a diameter of 320 fim, corresponding 
to 2.3 arcsec on the sky. A typical raw CCD t frame obtained with OPTOPUS 
with a scale of 9 A/pixel (resolution of 21 A) is shown in Figure 1. 

Figure 1 

There is a fundamental difference between long-slit and multifiber spectro- 
graphs: in the former, the sky spectrum is sampled on both sides of the target; 
in the latter some fibers are dedicated to the various targets, while the sky is 
sampled by a suitable number of different fibers, each with its own transmittance. 
In a multifiber spectrograph, in addition to the standard operations applied for the 
reductions of long-slit data, to evaluate the "true sky" it is necessary to determine 
the transmittance of each fiber and the infiuence of each fiber on the adjacent ones 
(cross-talk). 

t In the present observations a CCD ESO TEK 571 x 520, RON=8.8e-, 
GAIN=1.7e~/^i:>t/, pixels of 27//m was used. 
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Together with the "astronomical" frames, a few sets of cahbration frames 
are obtained: bias and dark exposures, dome flat-fields, internal-lamp flat-fields 
and wavelength calibration exposures. After the bias and dark correction of the 
frames, individual spectra are extracted integrating the flux along the columns 
of each channel * (7 pixels wide) centered by gaussian fitting of the transversal 
profile. For each frame this produces as many mono-dimensional spectra as the 
fibers in use. 

2.2) the "standard" method 

Following a "standard" approach, in order to determine the relative trans- 
mittance of the fibers we have used the fluxes recorded in the flat-field spectra, 
obtained with diffuse light on the dome (the results obtained with an internal lamp 
are deflnitely worse). 

For each sky spectrum we have calculated: 



- the symbol i indicates the channel (it identifles the spectrum and the flber too); 

- the symbols < S KY spectrurrii > and < FFi > are the values of the total flux 
into the channel between two given pixels yi and j/2 along the dispersion axis (e.g. 
50 - 450) for the sky and the flat-fleld spectra respectively; 

- the symbol < FF > is the average flirx of the observed flat-flelds < FFi > ; 

- < SKYi > represents the sky flux estimated from the flber. 

We can average over all the < SKYi >, obtaining < SKY > and compute: 



* The term channel in the following will indicate the sub-image of the CCD 
frame containing the signal coming from a given flber. 



< SKYi >=< SKY spectrurrii > • 



< FF > 

< FFi > 



where: 




< SKY, > 



- 1 



< SKY > 
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We can adopt the root mean square of Esky as an estimate of the accuracy 
of the sky subtraction. Such a procedure is imphcitly based on a number of 
hypotheses: 

- the shape of the spectral transmittance T(A) of the fibers is the same for aU the 
fibers: 

T,(A)/T,(A)^/(A) 

for each fiber i and j; 

- the influence (cross-talk) of each fiber on the adjacent ones is negligible and there 
is no stray light into the spectrograph; 

- the variations of the sky due to blue galaxies and airglow are negligible. 

The results obtained for some frames are shown in Table 1. Two frames for 
each field were obtained sequentially. The number of fibers dedicated to the sky 
varies between 14 and 36. 

Table 1 

As listed in Table 1 the accuracy of the sky subtraction is greatly variable 
(2 — 8%) from case to case, showing that the fiat-field based procedures provide 
typically poor and in any case unreliable results. 

2.3) the "5577 method" - an improved standard method 

In the spectral range accessible in our observations there are a number of sky 
emission lines, the most intense of which is the oxygen line [01] at 5577.4 A. 

The [01] line has the necessary qualifications to be used as an estimator of 
the fiber transmittance. In fact: 

1) its flux is sufficient to reduce the shot noise always below 0.6% and typically 
at 0.45% (with a telescope of 3.6 meters, a dispersion of 9 A/pixel and 40 
minutes of exposure time); 

2) at the faint magnitudes of interest in the present analysis, it is easily 
measurable both in the object and sky spectra; 
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3) it is subject to the same variations of the instrumental transmittance as the 
object flux (e.g. differential fiber flexures during the exposure, that however 
should not cause variations in the transmittance in a properly constructed 
system) . 

On the other hand there are also a few drawbacks: 

4) it represents the fiber transmittance only at 5577 A; 

5) it occupies only a small area on the CCD and consequently a "bad" pixel (a 
defective pixel, a cosmic-ray hit, etc.) can lead to a loss of accuracy; 

6) its availability is not guaranteed for all the configurations of the spectrograph. 

A subtle distinction between the "5577" and the "flat-field" methods has to 
be noted: while the variance of the latter is due only to errors in the calibration 
of the fiber transmittances, the variance of the former is due also to spatial sky 
variations. In this way the "5577" method potentially provides a more accurate 
sky subtraction. 

To estimate the accuracy of the sky subtraction we simply replace the flat- 
field flux < FFi > with the [01] flux < F5577i >, evaluated by fltting the line 
proflle along the dispersion direction with a suitable function (a gaussian turns 
out to be good approximation), to limit the influence of possible bad pixels. The 
results are shown in Table 2. 

Table 2 

The accuracy of the sky subtraction improves (to 2.1-3.0%) except for the 
frame C (shown in Figure 1) where an "anomalous" flber is present (see below); 
if this flber is rejected, the accuracy becomes 3.2% for this frame too. Above all 
the reliability of the results is greatly improved. 

3. ANALYSIS OF THE RESULTS 

We have investigated various hypotheses to justify the poor results provided 
by the "flat-fleld" method: 
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1) diffusion and reflections into the spectrograph; 

2) cross-talk among fibers; 

3) change from fiber to fiber of the shape of the spectral transmittance T(A); 

4) saturation and/or non linearity effects of the CCD; 

5) unknown dimensional dependent factors: possible correlations of the Esky 
with the number of the fiber i, with the position of the fibers on the focal 
plane of the telescope (right ascension, declination, distance from the center 
of the field, "azimuth" on the OPTOPUS template); 

None of the effects listed in the items 3, 4, 5 turned out to be sufficient to 
explain the improvement between the "fiat-field" and the "5577" method (for 
example Table 3 shows the effect of the variation of the transmittance as a function 
of the wavelength range estimated with dome fiat-fields; effects 4 and 5 have been 
measured to be even smaller). 

Table 3 



3.1 The Sky Concentration 

Let us address the effects of internal diffusion and refiections: we have grouped 
these phenomena under the name Sky Concentration or, in the following, SC. To 
extract the SC it is useful to study particular frames denominated mono- fiber flats 
in which only one fiber has been exposed to the dome or internal-lamp light. By 
two dimensional fitting with suitable functions it is possible to evaluate the SC 
contribution. Comparing a number of mono-fiber flats obtained exposing different 
flbers, it has been possible to deduce the following properties of the SC: 

1) it is independent of the particular flber exposed; 

2) it has rotational symmetry; 

3) its contribution is proportional (2.2%) to the total flux in the frame. 
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Figure 2 



After fitting the SC (Figure 2) with a two dimensional gaussian of ax = 
CTy = 95 ± 5 pixels the bias-subtracted astronomical frame have been corrected by 
performing: 

/ = O - NSC- < Ofiux > 

where : 

- /, O, NSC represent respectively, the corrected frame, the original frame and the 
flux normalized SC; 

- < Opiux > is the total flux of the original frame. 

The SC isn't responsible for the poor accuracy of the sky subtraction with 
the flat-fleld method. A flrst order computation applied on the {Esky)i expression 
demonstrates that, on frames on which the flbers are exposed at the same level, 
the SC has no effect. In fact, we can write: 



<SKYspectrumi>-<NSCi>-Fs 

(p N ^ <FFi>-<NSCi>-FF . 

{■t^SkyJi 1 <SKYspectrumj>-<NSCj>-Fs 



N ' 1^3=1 <FFj>-<NSCj>-FF 

where the new symbols < NSCi >,Fs,Fp represent the normalized Sky 
Concentration flux present in the channel, the total flux of the astronomical frame 
and the total flux of the flat-fleld frame, respectively. The sum is performed 
over all the N channels exposed on the sky. 

Assuming: 

- Fs 

< SKY spectrurrii >~ S ~ — 

Tp 

< FFi >~< FF >~ 

M 

where M is the total number of the flbers. Then: 
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The possible improvement is negligible and the application of the flat-field 
method to the SC corrected frames confirms this deduction. On the contrary, 
applying the "5577" procedure to the SC corrected frames, the sky subtraction 
accuracy improves to 1.6 — 2.4% (Table 3). 

3.2 The cross-talk influence 

After removing the SC from the mono-fiber flats, we have analyzed the 
extended broad tails in the transversal proflle of the mono-flber flats (Figure 
3). We have called transversal proflle, or briefly TP, a trace along a direction 
perpendicular to the direction of the dispersion (in our case approximately the 
same as the y axis). 

Figure 3 

A flrst-order argument as in the case of the SC, shows that the improvement 
in accuracy on Esky is negligible also in this case for frames with flbers exposed 
at roughly the same level. 

We concluded therefore that the source of unreliability of the "flat-fleld" 
method has to be ascribed to telecentricity effects and/or non- uniformity of the 
flat-fleld light source. 

4. THE INFLUENCE OF EACH FIBER ON THE ADJACENT 
ONES 

What are the effects of the cross-talk (in the following Point Spread Function 
Background or PSFB) on the "5577" method ? 

The broadness of the TP obviously indicates that also the PSF of the 
instrument has long tails at low flux levels. In principle it is possible to deconvolve 
the frame with a given PSF but problems arise (Brault and White, 1971), due 
mainly to the type of sampling and to the boundary conditions required (the 
signal doesn't go to zero within the boundaries of the CCD). We have preferred to 
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produce a simplified method by which we can subtract the cross-talk effect among 
the fibers. We assume that the PSF is a sum of functions of which only the first 
component (the main, PSFq) isn't responsible for the cross-talk. We can write: 

PSF(r) = PSFo(r) + PSFi(r) 

The PSF is assumed with a rotational symmetry. In fact all the emission lines 
of the helium lamp can be fitted by a two dimensional gaussian with the same 
standard deviation [ax = (Ty = 1.24 ± 0.05 pixels with the dispersion due mainly 
to differences among fibers), independent of the wavelength and the position of 
the channel on the CCD frame. In this way we can make use of the classical 
convolution theorem (Andrews and Hunt, 1977). The observed TP is the result of 
the convolution of the spectrum with the PSF: 

TP{x) = Spectrum PSF\\y^con8tant = TPo{x) + TP^{x) 

where: TPo{x),TPi{x) are due to PSFQ{r) and PSFi{r) respectively. It is useful 
to see that only the study of PSFi is important because only this part of the PSF 
produces the background PSFB. We find that the TPi can be described by a sum 
of gaussian functions: 

TP^{x) = J2Hk-e~^ 

k 

where: 

Hk, (7k are the amplitude and the standard deviations of the gaussian functions. 

If we hypothesize a PSFi as sum of gaussian functions, also TPi is a sum of 
gaussian functions (as it is observed) with the same standard deviations * : 

PSFi{r)^^Ak-e ''"^^ 

k 

According to the convolution theorem, the PSFB expression in the pixel space 

is: 

* The conservation of the function type doesn't hold for every function. 
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l-Xp + 1/2 /■j/p + 1/2 i-oo 

F{xp,yp)= dt du PSFi{t,u-y)J{y)- dy 

Jxp-l/2 Jyp-1/2 J -oo 

where: 

- J{y) is the flux distribution of the "true", ideal spectrum; 

- xp,yp are the coordinates of the pixel P; 

- F{xp,yp) represents the observed flux at the pixel P; 

which becomes with a first order development (for ak ^ I) and assuming the 
mono-fiber fiat positioned a,t x — 0: 

F{xp,yp) = 2^Ak-e ^ 'l^jiyi) -e ^--k' 

k i 

where j{yi) is the observed fiux of the spectrum at the pixel coordinate j/j; 

For a given yp = constant we obtain a particular TPi[xp). The coefficients 
Ak , are related by: 



In practice only 3 gaussian functions are necessary. In this way from the 
observed TP at a given yp it is possible to estimate by gaussian fitting the 
and the Hf. parameters (note that the Hj, depend on yp) from which the A^ 
(independent from yp) are obtained (Table 4). Then it is possible to estimate the 
overall background due to the PSFi by convolving an artificial frame (in which the 
fiux contained in each channel has been concentrated in one column corresponding 
to the X position of the center of the channel) , representing a rough approximation 
of a deconvolved image, with the PSFi (the main responsible for the cross-talk). 
This background PSFB (Figure 4) can now be subtracted from the SC corrected 
frame: 



II = 1- PSFB 
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where : 

- //, / represent respectively, the final frame and the SC corrected frame. 

Figure 4 

The improvement of the accuracy (from 1.6 — 2.4% to 1.2 — 1.6%, see Table 
2) depends strongly on the distribution of fluxes among the fibers of the original 
frame. The "anomalous" fiber in the frame C turns out not to be "anomalous" 
any longer. As shown in Figures 1 and 5 this fiber lies immediately to the right of 
a strongly exposed fiber, whose cross-talk was responsible for the "anomaly" . 

Figure 5 

Unfortunately the lack of information about the light flux outside the CCD 
frame causes a PSFB correction not entirely satisfactory, expecially in the red part 
of the spectra. The variations of the sky subtraction accuracy as a function af the 
wavelength range are shown in Table 5. 

Table 5 



5. THE EFFECTS OF IMAGE DISTORTIONS ON THE SKY 

SAMPLING 

With the "5577" method some defects of the "standard" approach have been 
eliminated. However, the subtraction of the sky from the object spectrum leaves 
embarrassingly large oscillations in the observed flux in correspondence to strong 
emission lines and bands. This problem strongly affects the red part of many 
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spectra making them in practice unusable beyond 7000 A. The usual spectra 
extraction procedure consists in the following steps: 

- spectra extraction, optimizing the Signal to Noise ratio, e.g. the procedure 
EXTRA/LONG in the MIDAS package (Banse et al., 1983); 

- calibration in wavelength of the spectra (sky and object spectra), rebinning 
with a constant wavelength step; 

- mean-sky evaluation; 

- sky subtraction; 

- absolute flux calibration. 

The second step partially removes the image distortions but introduces 
secondary efi^ects: not only the rebinning deteriorates the signal to noise ratio 
but also small errors in the wavelength calibration of each fiber give origin to 
large errors in the subtraction of the narrow sky emission lines. To overcome this 
problem we have developed a procedure which tries to map the image distortions 
in the pixel space. Using the Seidel distortion theory (Jenkins and White, 1957): 

p = Pq + E- pI + o(5) 

where: 

- p is the distance of a point from the optical axis on the object plane; 

- po is the distance of a point from the optical axis on the image plane; 

- £■ is the distortion coefficient of the optical system; 

- o(5) represents fifth order terms in po which are neglected in the third order 
theory. 

In an undistorted coordinate system the expression for a spectrum is: 

X = constant 

Its expression, in a first-order development, in the distorted coordinate system 
(the observed one) is: 
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x = ai + I3i-y + ^i-y'^ 



where: 

di = Qi + a — mh 
P^=m- 2hQiE 
li = QiE 

where: Qi are the x positions of the spectra at y = 0, m represents the tangent of 
the angle describing the misahgnment of the y columns of the CCD with respect 
to dispersion direction, a and h are the x and y coordinates of the optical axis on 
the CCD frame. 

By measuring ctj, 7j for each spectrum, we can solve the system and 
obtain the parameters a, 6, m, E (Table 6). 

Table 6 

The observed x and y coordinates of a given pixel are mapped on the corrected 
Fby: 

Y ^y + mx- E ■{y-h)-[{x-af + {y- hf] 

where mx represents the effect of the rotation (typically < 0.1 pixel, sec Table 6) 
and E ■ [y — h) ■ [{x — aY + {y — 6)^] the effect of distortion (it can easily exceed 1 
pixel). Only the Y coordinate is actually relevant because it represents the position 
along the dispersion direction of the resulting spectrum. After this correction the 
extracted spectra (monodimensional) show different dispersion scales and zero 
points, which can be estimated using the wavelength calibration exposures. The 
behaviour of the dispersion scales is shown in Figure 6. 

Figure 6 
14 



To bring all the spectra to the same wavelength scale and zero point it is 
necessary to introduce a further correction: 



Yci = coefi - Y + qi 

where: 

Yci indicates the corrected coordinate for the spectrum i ; 

coefi , Qi are the dispersion scales and the zeropoint shift values of the i spectrum, 
estimated by comparing the positions of He-Ar emission lines with a reference 
spectrum. From Figure 6 it is possible to deduce that the effects due to the 
different scales imply displacements of the order of half a pixel. 

To determine the residual shifts, essentially due to differences between the 
optical path of astronomical observations and He-Ar calibrations (Figure 7, bottom 
panel), we can use a few emission sky lines (e.g. the [01] 5577 A, 6300 A) and 
carry out the final correction: 

Yci = Yci + Ti 

where: 

Yei are the final corrected coordinates for the spectrum z*'*; 

Ti is the shift of the i^^ spectrum. 

After all these corrections the sampling of each spectrum is difi'erent from 
the other ones. The union of all the sky normalized spectra (the SKYi spectra) 
allows the building of an "oversampled" global sky spectrum. By interpolation 
and rebinning of this "oversampled" spectrum we obtain an estimate of the sky in 
the reference system of wavelengths of each object spectrum. 

Figure 7 

The advantages of this sky subtraction technique are: 
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- the drastic reduction of the oscillations in the sky subtracted spectra in 
correspondence of narrow emission features 

- preservation of the true signal to noise ratio and of the independence of the 
flux of each pixel, allowing also a correct estimate of the S/N 

Typical results are shown in Figure 8 and 9. 

Figure 8 and 9 

The procedures described in sections 4 and 5 are illustrated in the flow chart 
of Figure 10, representing a set of procedures developed in the framework of the 
MIDAS image processing system (Banse et al, 1983). 

Figure 10 

6. CONCLUSIONS 

The poor and unreliable sky subtraction often observed when operating a 
multifiber spectrograph depends strongly on the method used to estimate the 
transmittance of the fibers. A method based on the observed flux of the 5577 [01] 
line has been shown to allow a considerable increase in the accuracy of the sky 
subtraction. 

Further improvements are obtained when the effects of the internal scattered 
light and the influence of each fiber on the adjacent ones are determined and 
corrected. 

For the OPTOPUS spectrograph, this allows to reach a sky subtraction 
accuracy of 1.2-1.6%. 

To obtain the maximum performance from the instrument it is necessary to 
develop optimized extraction and subtraction procedures that take into account 
also the effects of the image distortions on the sky sampling. 
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The overall improvement is important in view of the application of multifiber 
instruments on telescopes of 6-10 meters diameter (VLT, KECK, etc.). With low 
dispersion (~ 10 A/pix ) and reasonably long exposures (>~ 1 hour) the main 
limitation to the S/N already for telescopes of the 4m class is due to the sky 
subtraction, if its accuracy is limited to the 5-8% range. In such a situation it 
would not be useful in terms of limiting flux to apply multi-fiber spectrographs 
on bigger telescopes. The techniques described in this paper demonstrate that 
this limit can be significantly pushed down, making worthwhile the construction 
of multi- fiber spectrographs for very large telescopes. 

ACKNOWLEDGEMENTS: We thank P. Andreani, G.Avila and S. 
D'Odorico for enlightening discussions and suggestions. 



17 



REFERENCES 

Andrews, H.C. and Hunt, B.R., 1977 Digital Image Restoration, Prentice-Hall 
Signal Series. 

Avila, G. and D'Odorico, S., 1993 Fiber Optic in Astronomy, ASPC. 37, 90. 

Banse, K., Crane, P., Ounnas, C. & Ponz, D., 1983. In: Proc. of DECUS, Zurich, 
p. 87. 

Brault, J.W., White, O.R., 1971, A&A, 13, 169. 

Elston, R. and Barden, S., 1989, NOAO Newsletter, 9/89, p. 21. 

Jenkins, F.A., White, H.E., 1957, Fundamentals of Optics, ed. McGraw-Hill, New 
York, pg. 132 

Mignoh, M. and Cuby, J.G., 1994, Proc. of SPIE conference 2198, in press. 

Schwarz, H.E. and Melnick, J., 1993 ESO Users Manual. 

Wynne, Ch.G., 1993, MNRAS, 260, 307. 

Wyse, R.F.G. and Gilmore, G., 1992, MNRAS, 257, 1. 



18 



FIGURE CAPTIONS 



Figure 1: A typical raw CCD frame obtained with OPTOPUS (in this case 
frame C). Notice the 22"''^ fiber from the left, strongly exposed at the ambient 
light. The "anomalous" fiber discussed in the text lies immediately next to the 
right. 

Figure 2: Representation of the Sky Concentration (in units of 10~^). 
Figure 3: The mean Transversal Profile (see text). 

Figure 4: The Point Spread Function Background (PSFB) for the frames D 
and A. Notice the effect of the one strongly exposed fiber in frame D. 

Figure 5: The calculated PSFB at row 300 of frame C. 

Figure 6: The variation of the dispersion coefficient versus the coordinate of 
each channel. 

Figure 7 : Effect of the residual shifts qi between the calibration and object 
frames (frame B). Top: all the sky spectra in the Yci system coordinate (see text). 
Middle: the same spectra in the Yci system coordinate. Bottom: the shift values 
versus the center coordinate of the channel along the x axis. 

Figure 8: Spectra of the same object extracted with the "standard" reduction 
(top) and with the procedure described in the text (bottom). 

Figure 9: A new quasar (continuous line) discovered in the field of SA94, 
compared with the sky flux (dashed line): only with the new procedure it has been 
possible to extract and identify the characteristic emission lines of this spectrum. 

Figure 10: Flow chart of a reduction session of OPTOPUS frames using 
the procedures produced within the framework of the MIDAS image processing 
system. 
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